Mass spectrum and Levy— Schrodinger relativistic equation 
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We introduce a modification in the relativistic equations in such a way that (1) the relativistic 
Schrodinger equations can always be based on an underlying Levy process, (2) several families of 
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I. INTRODUCTION AND NOTATIONS 

In this note we adopt the space-time relativis- 
tic approach of Feynman's propagators (for bosons 
and fermions) instead of the canonical Lagrangian- 
Hamiltonian quantized field theory. Indeed the former 
alternative is preferred to the latter for the developments 
of our basic ideas that exhibit the connection between the 
propagator of quantum mechanics and Levy's stochastic- 
ity. More precisely the relativistic Feynman propaga- 
tors are here linked with a dynamical theory based on a 
particular Levy stochastic process. This point, already 
mentioned in a previous paper [l[ , is here analyzed thor- 
oughly with the purpose of deducing its consequences for 
the case of fundamental fermions and bosons (quarks, 
leptons, gluons etc. . .) of the Standard Model (SM) char- 
acterized by the symmetry SU C (S) x SU L (2) x £7(1). To 
this end we now recall that a Levy process is a stochastic 
process X(t), t > on a probability space (fi, F, P) such 
that 



X(0) = 0, 



-q.o. 



• X(t) has independent and stationary increments: 
for each n and for very choice of < t\ < t<z < 
. . . < t n < +oo the increments X(tj+i) — X(tj) are 

independent and X(tj+i)— X(tj) = X(tj+\—X(tj); 

• X(t) is stochastically continuous: for every a > 
and for every s 

\imW (\X(t) - X(s)\ >o)=0. 

t— >s 



To simplify the notation we will consider in the following 
one-dimensional (the n-dimensional extension would not 



be a difficult task) Levy processes: it is well known 0- 
4] that all its laws are infinitely divisible, but we will 
be mainly interested in the non stable (and in particu- 
lar non Gaussian) case. In other words the character- 
istic function of the process At-increment is [f(u)] At ' T 
where tp is infinitely divisible, but not stable , and r is 
a time scale parameter. The transition probability den- 
sity p(2 1 1) = p(x2, t?.\x-i, t\) of a particle moving from the 
space-time point 1 to 2 then is 

1 r +co 

p(2\l) = — du[<p{u)] {t2 - tl)/T e- iu{x *- Xl) (1) 

27r J-oo 

In analogy with the non relativistic Wiener case we ob- 
tain for the motion of a free particle the Feynman prop- 
agator /C(2| 1) = K.(x%, t%\xi, ii) as 



1 r + °° 

£(2|1) = — / du\<p(u)] t(t2 - tl)/T e- tu{x2 - Xl) 
2tt J-oo 

and the corresponding wave function evolution is 

r + CO 

ip(x,t) = / dx' JC(x, t\x' , t')tp(x' , t'). 



(2) 



(3) 



From ([2]) and ((3]) we easily obtain [l[ 
idt-ifi = --r)(d x )ip 

T 

where rj = log(p and r/(d x ) is a pseudodifferential opera- 
tor with symbol rj(u) defined through the use of Fourier 
transforms 0, 043 ■ It plays the role of the generator of 
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1 A law ip is said to be infinitely divisible if for every n it exists a 
characteristic function ip n such that tp = i^™; on the other hand 
it is said to be stable when for every c > it is always possible 
to find a > and b € R. such that e lbu tp(au) = [ip(u)] c . Every 
stable law is also infinitely divisible. 



2 



the semigroup T t = e t11 ( 9x " r operating on the Banach 
space of the measurable, bounded functions 

It is very well known [2,01, on the other hand, that tp 
represents an infinitely divisible law if and only if rj{u) = 
\ogtp(u) satisfies the Levy-Khintchin formula 



rj(u) = i'yu 



8 2 n 2 

-T—+ / [e lux - 1 -iux /[_!,!] (x)] u(dx) 

(4) 

where 7, /? € M, I[_ 11 ](x) is the indicator of [—1, 1], and 
v(dx) is the Levy measure, namely a measure on K such 
that v({0}) = and 



(x 2 A 1) v(dx) < +00. 



In the case of a centered, symmetric law the equation ((4]) 
simplifies in 



r)(u) 



2, ,2 



p 2 u 



+ / (cosux — 1) v(dx) (5) 



and rj{u) becomes even and real. As a consequence the 
corresponding operator r](d x ) is self-adjoint and acts on 
propagators and wave functions according to the Levy- 
Schrodingcr integro-differential equation 



and for a plane wave 4*{ x ) 
metric Levy noise 



from ([5]) with a sym- 



E n <b(x) = — 



p 2 u 2 
2 

P 2 u 2 



+ / (e ±iuy - 1) v{dy) 
Jn 

+ / (cosuy — 1) v{dy) 



4>{x) 



which is satisfied when Eq = —arj{u)/T. Finally by tak- 
ing p = au for the momentum we obtain the relevant 
equation 



E 



a (p 
--77 1 - 
r \a 



(9) 



III. RELATIVISTIC QUANTUM MECHANICS 

Equation © connects the kinetic energy of a force- 
less particle to the logarithmic characteristic of a Levy 
process: if we take in particular the non stable law 
rj(u) = 1 — Vl + a 2 u 2 with the following identification 
of the parameters 



= K 



= mc , 



p = hu. 



id t ip(x,t) 



1 



r)(d x )ip(x,t) 



(6) 



1 



[ip(x + y 1 t)- %jj(x, t)} v{dy). 



The integral term accounts for the jumps in the trajecto- 
ries of the underlying stochastic process, while an action 
a with (3 2 = ar/m provides the usual differential term 
of the Schrodingcr equation. For (3 — a pure jump 
Levy-Schrodingcr equation is obtained 



id t ip(x,t) 



1 



[ip(x + y,t)-ip(x,t)} v{dy). (7) 



II. STATIONARY SOLUTIONS FOR THE FREE 
PARTICLE 

Equation ([6]) allows a simple stationary solution: if we 
consider 



to/3 2 



i/>(x,t) =e-' lEot/a (j)(x), 
we have then 



E <t>(x)=-?L<f>"(x)-" I [<f>(x+y)-<f>(x)}v(dy), (8) 



we are led to the formula 

Eq = —mc 2 r] (j^j — E — mc 2 = \J m 2 c 4 + p 2 c 2 — mc 2 

(10) 

which is the well-known relativistic kinetic energy for a 
particle of mass to. The Schrodinger equation of a rela- 
tivistic free-particle is easily obtained from (|10|) by rein- 
terpreting as usual E and p respectively as the operators 
ihdt and —ihd x : 



ihd t ip(x,t) = v/to 2 c 4 - h 2 c 2 d 2 ip(x,t), (11) 

but this comes also from © after having absorbed the 
mass energy term —mc 2 of (flU)) into a phase factor 
e lmc t / h . In three dimensions (TTTj) reads 



ihd t ip{x,t) = ^m 2 c i - h 2 c 2 V 2 ip(x,t) (12) 

It has been shown 0, H{ that the Levy process behind the 
equations (fTTj) and (fT2)) is a pure jump process [l|, Q with 
an absolutely continuous Levy measure v(dx) = W{x)dx 
and 



7r|x| 



V h 



(13) 



(K v are the modified Bessel functions @), that in three 
dimensions becomes 



W(x) 



1 



2a7r 2 | 



, \x\ 



2tm 2 \ 



K 2 —\x 



(14) 
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while from (JT)) the equation (fTTj) becomes equivalent to 

ihdtip(x,t) (15) 
2 f ip(x + y,t)-ip(x,t) (mc \ 

= ~ mc i 1 VT M J dy 

and in three dimensions 

ihd t ip(x,t) (16) 
2 /" ip(x + y,t) — tp(x, t) mc , /mc \ 3 



-me 



27r 2 |y| 



t^ 2 It |y| ' y 



From the equation (j!2j) by the well known standard pro- 
cedures [l3| one derives (always for the free particle) the 
Klein-Gordon and Dirac equations in three dimensions 
for the wave functions and spinors t/j, respectively 



□ 



n 2 



»7m 9M - -J-) i> 



i> = 0, 
0. 



(17) 
(18) 



The Klein-Gordon and Dirac propagators verify instead 
the inhomogencous equations (here h — c = 1) 



(□2 -■ 



2 )K, KG {2\\) = <5 4 (2|1) (19) 
m)JC D {2\\) = i5 i {2\\) (20) 



with <5 4 (2|1) = 6(t 2 - h)S 3 {x 2 - an). Let us finally re- 
mark that these relativistic quantum wave equations have 
been recently of particular interest [ll| also in the field 
of quantum optical phenomena and of quantum informa- 
tion. 



IV. INFINITE DIVISIBILITY-PRESERVING 
MODIFICATIONS 

The relativistic, time-like four-momentum p = (E/c,p) 
obeys the relation (here p 2 will represent the square mod- 
ulus of the tri-vector p) 



2 - L/ 2 2 2 ^ n 
p = —5 — p = m c > 



(21) 



so that the hamiltonian dependence of energy on momen- 
tum is 



E(p) = VmV + P 2 c 2 = mc 2 J 1 + • (22) 

V ?71 Z C Z 

This relativistic energy E containing a rest mass term 
mc 2 , the kinetic energy E Q in a dimensionless form be- 
comes 



(23) 



By taking now 



E P 

77= M= 

mr amc 



where a is a constant with the dimensions of a length, 77 
will be dimensionless while u will be the reciprocal of a 
length, and the equation (|23"1) becomes 



77(11) = 1 — vHoV 

namely - not surprisingly - it takes the three-dimensional 
form of the logarithmic characteristic giving rise to the 
relativistic quantum equations in the Section [IIII 

Our purpose consists now in proposing a modification 
of 77(11) that preserves the infinite divisibility of the law, 
and eventually produces changes in the forceless equa- 
tions of motion for a particle - both at the classical and 
at the quantum level - in comparison with the classi- 
cal and quantum motions given by the equations (|17[) 
and (fT5|) . To this end we modify (|2"2"|) in the following 
way 



E(p) 



mc 




P 



+ f 



P 



(24) 



where / is a - possibly small - dimensionless, smooth 
function of the relativistic scalar p 2 /m 2 c 2 . Of course this 
modification entails that p 2 no longer coincides with m 2 c 2 
since the standard relation (|21[) is now changed into 



P 



p2 

2 2 2, 2 2 c 

— p — m c +m c j 



p 



(25) 



As we will see in the following this also implies that the 
mass no longer is m: it will take instead one or more 
values depending on the choice of /. In fact it could ap- 
pear to be preposterous to introduce a function / of an 
argument which after all is a constant (albeit different 
from 1). However we will show that this will lend us the 
possibility of having both a mass spectrum, and a new 
wave equation when - in the next section - we will quan- 
tize our classical relations. Moreover it will be argued in 
the following that in this way the corresponding modified 
logarithmic characteristic 77 will remain infinitely divisi- 
ble: a feature that is instrumental for a connection to a 
suitable underlying Levy process. 

To see that we first remark that (|25[) defines the total 
particle energy E in an implicit form. To find it explicitly 
we rewrite (|25p in a dimensionless form as 



1>- 



P 



1 + / 



and then, by taking g(x) = x — f(x), we just observe 
that the former equation requires that x be solution of 
g(x) — 1, namely 



V 



P 



f 



P 



1. 



If then g represents one of the (possibly many) so- 
lutions of this equation, we could write 



4 



so that wc have 



E 2 

p 2 = — -p 2 = m 2 c 2 g- 1 (l) 



which can be interpreted as a simple mass re-scaling from 
m to one of the (possibly many) values M = m-v/ <7 _1 (1). 
The new hamiltonian then is 



E(p) = v / m 2 c i g~ 1 (l) +p 2 c 2 = Mc 2 \ 1 



M 2 c 2 



(26) 



and its kinetic part (by applying the same re-scaling also 
to the subtracted rest mass term) is 



Eo(p) = E(p) - mc 2 y/FHTj 
= Mc 2 



1 



P -Mc 2 . 



M 2 C 2 



Hence the main consequence of our modification consists 
of a re-scaling of the mass value (m — > M) at a purely 
classical level. This fact is apparently helpful because 
it is straightforward to see now that the new associated 
logarithmic characteristic 77 is again infinitely divisible, 
and hence still produces acceptable Levy processes. But 
there is more: since c/ _1 (l) can take several different real 
and positive values, by means of our modification (|24|) 
we have introduced an entire mass spectrum: indeed in 
the rest frame of the particle we have 



M = E cm /c 2 = my/g-^l) 



(27) 



V. QUANTUM EQUATIONS OF MOTION 

It is important to remark that while the equation (|24[) 
allows a peculiar transition to quantum mechanics (E — > 
ihdt, p — > — iTiV) if we interpret this energy formula as 
a new hamiltonian operator, namely it leads to 



ihdtip(x,t) = mc 2 



h 2 



rV 2 + / 



□ 



ip(x,t) 
(28) 

the equation (|26| gives instead the usual Klein-Gordon 
equation (jT2"j) with just a possibly re-scaled mass M = 
g~ l {\). In fact from (|28[) one obtains now a modi- 
fied Klein-Gordon equation for both the wave function 
ip and its corresponding propagator JCkg(^\^) (from here 
on h = c = 1) 



□2 - m 2 f 



1 



1 



□ 



■;?).- 







ICkg(2\1) 



(29) 



(30) 



6\2\1) 



and by standard methods [T(J the modified Dirac spinor 
equations 



□ 



ip = (31) 



i^-mJl + fl-sD 



£ (2|1) (32) 
= iS 4 {2\l) 



In the momentum space (with Fourier transforms in four 
dimensions) these equations become much simpler: more 
precisely we have 



£kg(p 2 ) = 

ic d (p 2 ) = 



1 



p2 _ m 2 _|_ j(jp. jm 2 )] _|_ jg 
1 



l^Pn — m-y'l + f(p 2 /m 2 ) + it 



We notice that /Crj(2|l) is in our case simply related to 
the /Ck*g(2|1) (like in the usual case) as 

JC D (2\1) = * (i? 2 + my/1 + f(a 2 /m 2 )) JC KG (2\1) 



VI. PHENOMENOLOGY: QUARK AND 
LEPTON MASSES 



The equations ([30]) and (|32|) generalize the well known 
propagator equations (|19|) and (|20[) which derive from 
QED and QCD at zero order (in absence of interaction 
terms). For future developments we recall that the La- 
grangian density of QCD is 2 , up to gauge fixing terms: 



where = d^A^ - d u A a ^ + gsfabc^A^,, and the inser- 
tion of interaction terms is done with the minimal inter- 
action by substituting the simple derivative with the 
covariant one where we have respectively for QED 
and QCD 

Dp = dfi - ieA^ 

(£>„)<. = 6,a, in-i/A';,. 

The Standard Model (SM) SU C (3) x SU L (2) x U(l) treats 
both strong and electro-weak interactions: within this 
scheme the modified rj(u) leads to new interesting conse- 
quences. We begin by considering the Feynman rules in 



Here g a is the QCD coupling constant, and / a j, c are the 
SU(3) color matrices and structure constants respectively, and 
A® the eight Yang-Mills gluon fields; i/>? are the Dirac 4-spinors 
associated with each quark field of color i and flavor q. 
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perturbation theory in presence of the modified zero or- 
der propagator for both spin i (quarks and leptons) and 
spin 1 (gluons, vector weak interacting Bosons). The 
amplitude A for a fermion that propagates from ver- 
tex X to vertex Y if expanded looks as follows: A = 
A(°) + A^ + A^ + ... The lowest order is 



AW = Y- 



-X. 



j^Pfj, - my/1 + f(p 2 /m 2 ) + ie 

It is then possible that the Fermion emits and reabsorbs 
a virtual vector boson from X to Y: 



A^ = 4ng 2 Y / d 4 k 



r 



YPp - my/1 + f(p 2 /m 2 ) (p - k) 2 
1 



k v ^ v — my/1 + /(fe 2 /m 2 ) + ie 



In 



X 



■yPpp - my/1 + f(p 2 /m 2 ) 

We choose now fix) in such a way that it makes finite 
the integral 



d 4 fc 



jPkp - my/1 + f (k 2 /m 2 ) + ie (p - k)'- 



— 7T^> 
(33) 

One may notice that fix) behaves as a smooth cut-off 
in a procedure of Regularization at each order in QCD 
(and QED). The integral C is an invariant of the form 
C = Aip 2 )$ — Bip 2 ) and its integrand is also present as a 
factor in higher order terms, thus producing convergence. 
In a similar way one expects that the representation of 
the complete fermionic propagator, as well as its zero- 
order, is made up of two additive terms, in momentum 
space, each of them exhibiting simple analyticity prop- 
erties except for a limited number of poles and branch 
points 3 [l4|: for the zero order we have 



+ my/l + fjp 2 /m 2 ) 
! — m 2 [l + fip 2 jm 2 )\ 



$ — my/l + fip 2 /m 2 ) 



Let us now reconsider the equation ([55)) . The simplest 
expression for f(x) compatible with C finite is a polyno- 
mial of third degree in x: 



fix) = \ix + \ 2 x 2 + X 3 x 3 



(34) 



and this is suggestively connected with the possibility 
of having a mass spectrum. Indeed, as stated in the 
Section IIV1 the spectrum is produced by the multiple 
solutions of the equation <?(x) = x — fix) = 1, and 
we achieve three values that, under proper conditions, 
might correspond to three masses. If we consider the 



3 Likewise one can reason for the complete and zero order propa- 
gator of basic bosons (gluons, W^, Z°, Higgs). 



m d 


m 3 


m b 


3 x 10-* 


70 x 10-* 


4.13 


7 x Kr a 


120 x 10~ a 


4.27 


Ai 


A 2 


A 3 


-1.84 x KT 3 


1.84 x 10~ a 


-9.69 x 10~ ilJ 


-3.41 x 10" a 


3.41 x 10" a 


-9.14 x 10~ a 



TABLE I: Estimated values of the A's in 
charge — |. Masses are in Gev/c 2 . 



for quarks with 



m u 


m c 


m t 


1.5 x 10~' 5 


1.16 


171.2 


3.0 x 10- 3 


1.34 


174.0 


Ai 


A 2 


A 3 


-1.67 x 10-" 


1.67 x 10~° 


-1.28 x 10~ ib 


-5.01 x 10"° 


5.01 x 10~ b 


-1.49 x 10- Lb 



TABLE II: Estimated values of the A's in ([34j) for quarks with 
charge |. Masses are in Gev/c 2 . 



m e 




m T 


5.11 x 10~ 4 


105.6 x 10 _ 3 | 


1.77 


Ai 


A 2 


As 


-2.35 x 10 _b 


2.35 x 10 _b | 


-1.95 x 1Q- 12 



TABLE III: Estimated values of the A's in (01} for leptons 
with charge —1. Masses are in Gev/c 2 . 



three (real and positive) zeros x\,X2 and x% of the poly- 
nomial gix) — 1 = x — fix) — 1 we easily find the following 
simple algebraic relations with the A's: 



Ai 
A 2 



1 



1 

X2 



X3 
1 



A, 



XiX 2 XiX 3 x 2 x 3 



XlX 2 X 3 



The connections with the possible experimental physical 
masses are Mi = my/xi, M 2 = my/x^, M 3 = 
If the three poles in the free (zero order) propagator are 
real and positive (with proper residues) , with appropriate 
values of the A's, they allow the interpretation of physical 
basic masses of fermions (quark or leptons) belonging 
to the three different families of the Standard Model. 
To be more specific we get two different propagators for 
quarks, one with charge — | (d,s,b quarks) and another 
with charge +| (ti, c, t quarks). Similarly for charged 
leptons (charge —1 and spin ^) we get one propagator. 

In the Tables HI [TT| and IIIII we give a few examples of 
numerical values of the A parameters in equation (|34l) for 
quark and lepton masses (in Gev/c 2 ) taken from the Par- 
ticle Data Group [l3| • For the quarks with charge — -| and 
different mass estimates we obtain the results listed in 
the Table HI while for quarks with charge | we get the re- 
sults of the Table HH and finally for charged leptons with 
charge — 1 we have the results of Table IIIII According to 
our model a significant contribution to the fermion mass 
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spectrum derives from the poles of the zero order prop- 
agator, whereas the role of the interaction terms might 
become complementary and can be estimated in (renor- 
malized) perturbation theory. 

Our A's also produce the regularization at each order 
in QCD (and QED) perturbation theory. Furthermore 
we know that after regularization the approximate rep- 
resentation of the propagator tends to a finite limit (ex- 
act propagator) due to the renormalization mechanism. 
More precisely in these field theories the calculated renor- 
malized physical quantities are supposed to become inde- 
pendent of the cut-off. The latter must disappear in the 
transition from regularization to renormalization. Conse- 
quently the A's remain in a finite fixed number and tend 
to definite (real) values as the cut-off cancels out. Within 
this scenario the poles representing physical masses (3 in 
our case) remain stable (even if shifted partly with re- 
spect to those computed approximately by us), because 
of to the assumed analytical properties of the renormal- 
ized propagator [l2|, [H| • 

VII. CONCLUSIONS 

We have proposed a modification of the classical rel- 
ativistic hamiltonian that allows the presence of several 
masses without changing its basic structure. This mod- 
ification does not affect the infinite divisibility of the 



laws that are at the basis of the correspondence between 
stochastic processes and Levy-quantum mechanics equa- 
tions. However we discovered that the mentioned modi- 
fication suggests a reformulation of the relativistic equa- 
tions for wave functions and propagators in such a way 
that a suitable choice of the background noise produces 
a convergence in the perturbative contributions. To this 
purpose we remarked that a modification - with respect 
to the one given by equation (fT0|) - of the logarithmic 
characteristic r](u) by the insertion of the cut-off f(x) al- 
lows to proceed to Regularization first, and then Renor- 
malization of the two-point function of QCD. There are 
three parameters in our phenomenological f{x) which is 
a third degree polynomial; the latter appears as the sim- 
plest choice that produces convergence in the integrals 
representing high order contributions to the fcrmion and 
boson propagators. Such parameters create three dif- 
ferent poles in the zero-order propagators and allow the 
interpretation of a physical system with three different 
masses under precise constraints on f{x). The masses 
might be related to the three families of the Standard 
Model. We like to point out that from the analytic- 
ity properties of the renormalized propagators, the men- 
tioned poles tend to stabilize in the limit toward the com- 
plete solution (even if shifted with respect to the zero 
order ones) whereas the cut-off is expected to disappear 
because of the regular renormalizable QCD theory fl2l |. 
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